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ABSTRACT 

We report the discovery with the 100 m Green Bank Telescope of 65 ms radio pulsations from the X-ray pulsar 
J0205+6449 at the center of supernova remnant 3C58, making this possibly the youngest radio pulsar known. 
From our observations at frequencies of 820 and 1375 MHz, the free electron column density to PSR J0205+6449 
is found to be 140.7 ± 0.3 cm -3 pc. The barycentric pulsar period P and P determined from a phase-coherent 
timing solution are consistent with the values previously measured from X-ray observations. The averaged radio 
profile of PSR J0205+6449 consists of one sharp pulse of width « 3 ms « 0.05P. The pulsar is an exceedingly 
weak radio source, with pulse- averaged flux density in the 1400MHz band of ~ 45/iJy and a spectral index of 
~ -2.1. Its radio luminosity of ~ 0.5 mJy kpc 2 at 1400MHz is lower than that of ~ 99% of known pulsars and is 
the lowest among known young pulsars. 

Subject headings: ISM: individual (3C58) — pulsars: individual (PSR J0205+6449) — supernova remnants 


1. INTRODUCTION 

The “Crab-like” supernova remnant (SNR) G1 30.7+3.1 
(3C58) has long been suspected of containing a pulsar at its 
center (e.g., Becker, Helfand, & Szymkowiak 1982). This SNR 
is of particular interest because it is thought to be the remnant 
of the explosion recorded in 1181 AD (SN 1181), and hence to 
have a known age of 820 yr. After more than 20 years of search- 
ing in the X-ray and radio bands, the pulsar J0205+6449 with 
period P = 65.68 ms was recently discovered at the SNR center 
with the Chandra X-ray Observatory (Murray et al. 2002). 

The historical record associating the pulsar-powered syn- 
chrotron nebula 3C58 with SN 1181 is compelling. Neverthe- 
less, 3C58 has some properties suggesting a larger age (see, 
e.g., Slane, Helfand, & Murray 2002, and references therein). 
Also, the nebular radio emission of 3C58 is ~ 10 x lower than 
the comparably aged Crab Nebula’s, while the X-ray emission 
is ~ 1000 x weaker. Because the synchrotron lifetime of X- 
ray-emitting electrons is short (unlike of those generating radio 
emission), it is surprising that the ratio of the X-ray luminosities 
from the 2 SNRs does not match the relative present-day spin- 
down luminosities of the pulsars: £crab/ £ jo 205+6449 w 15. Since 
E oc P/P 3 , the power injected into the nebula as a function of 
time depends strongly on P(/), which therefore is an important 
quantity to constrain; the subsequent nebular evolution is gov- 
erned in addition by the pressure of the ambient medium. The 
characteristic age r c = P/2P = 5400 yr of J0205+6449 can be 
reconciled with the 1181 AD explosion by appealing to a large 
initial spin period of ~ 60ms (Murray et al. 2002). However 
we do not know the form of the spin evolution of the pulsar: 
measuring its “braking index” /?, where the angular frequency 
of rotation evolves as ft ex -ft /l , would constrain P(t) and help 


to understand the energetics of 3C58. 

The study of J0205+6449 at radio wavelengths could con- 
tribute significantly towards the determination of the braking 
index: precise long-term monitoring of the pulsar rotation, es- 
pecially important for young pulsars which often experience pe- 
riod glitches, may prove easier from the ground than from satel- 
lites. In addition, the absolute phase alignment of radio and X- 
ray pulses should provide information regarding the emission 
mechanism. A radio detection would also yield the dispersion 
measure (DM), giving an independent distance constraint or, 
alternatively, the average free electron density along the line of 
sight. Finally, measuring the radio luminosity and spectrum of 
very young pulsars is essential for constraining the population 
of these objects in the Galaxy. Motivated by these factors we 
made deep searches for radio pulsations from PSR J0205+6449, 
culminating in their discovery with the new Green Bank Tele- 
scope as reported here. 

2. OBSERVATIONS 

The SNR 3C58 was most recently searched for a radio pulsar 
by Lorimer, Lyne & Camilo (1998), with a resulting flux den- 
sity limit for the pulsar of $600 <1-1 mJy at a search frequency 
of v = 600 MHz. This corresponds to Shoo < 0.3 mJy for a 
median pulsar spectral index of a = -1 .6 (Lorimer et al. 1995), 
where S u oc v a . At a distance d = 3.2 kpc (Roberts et al. 1993) 
the implied luminosity limit Lu oo = Smoo ^ 2 < 3mJykpc 2 is 
somewhat larger than the radio luminosity of at least one young 
pulsar (Camilo et al. 2002), and the discovery of X-ray pulsa- 
tions from PSR J0205+6449 at the center of 3C58 by Murray et 
al. (2002) led us to attempt more sensitive searches. 

On 2001 October 8 we used the 100 m Effelsberg telescope to 
conduct a search of PSR J0205+6449. We observed the pulsar 
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position for 12.5 hr at a center frequency of 1410MHz. Sig- 
nals were filtered in 8 x 4 MHz-wide channels for each of 2 
polarizations, after which they were detected and digitized ev- 
ery 0.8 ms. Samples from orthogonal polarizations were added 
in hardware and the 8 total power time series were written to 
disk for offline analysis. We also recorded short data sets on 
strong pulsars for calibration purposes. We analyzed the data 
in standard fashion (see below for description of the search at 
Green Bank), de-dispersing the time series for a number of trial 
DMs in the range 0-500 cm -3 pc, with no significant candidates 
found at the expected period. From these observations we in- 
ferred 5*1400 < 0. 1 mJy for a pulsar duty-cycle of ~ 0.05F. 

On 2002 February 22-23 we used 18 hr at the Green Bank 
Telescope to search for radio emission from PSR J0205+6449. 
We divided the time nearly equally between an observation cen- 
tered at 1375 MHz with a Gregorian focus receiver, followed by 
one at 820 MHz with a prime focus receiver. The telescope gain 
and system temperature, including a large contribution from the 
SNR, are given in Table 1. We also performed short calibra- 
tion observations of strong pulsars at both frequencies. The re- 
ceivers, sampling 2 orthogonal polarizations, had band-limiting 
filters installed of width 150 and 80 MHz centered on 1375 and 
820 MHz respectively. The radio frequency (RF) signals were 
converted to an intermediate frequency (IF) and transmitted via 
optical fibers to the control room, where they were passed to 
the Berkeley-Caltech Pulsar Machine (BCPM). 

The BCPM is one of several signal processors of similar de- 
sign (Backer et al. 1997). The 2 input IF signals are distributed 
to an array of 2 x 6 mixer/filter modules. These accept sig- 
nals from 6 local oscillator modules for baseband mixing and 
low-pass filtering with the bandwidth set by subsequent signal 
processing requirements. The signals are sent to 2 x 6 digital 
filter boards which 4-bit sample and form detected powers in 
16 independent and adjacent frequency channels. The BCPM 
is thus an analog/digital filter bank with 2 x 96 channels. Time 
decimation of the power samples and optional summing of the 
2 IFs (which we used) is accomplished in combiner boards. The 
combiner board output is reduced to 4-bit power deviation by 
mean removal and scaling. The power vector(s) are passed to 
the host Sun workstation through an EDT 9 DMA card. The Sun 
attaches a header and places data on a disk for offline process- 
ing as well as monitoring data quality 10 . For our 2 observations 
we used the data-taking parameters listed in Table 1. 

3. DATA ANALYSIS AND RESULTS 

Although we knew the pulsar period a priori, we did not 
know its DM and chose to use standard search algorithms to 
analyze the data. Also, despite the presence of strong RF inter- 
ference (RFI) visible in a spectrum analyzer during data collec- 
tion (most notably as persistent radar in the 1400MHz band and 
sporadic but very strong RFI of unknown origin in the 800 MHz 
band) in a first pass we masked no portions of the fluctuation 
spectra. For both observations data from 2 25 decimated time 
samples (see Table 1) were de-dispersed at 204 trial DMs in 
the range 0-203 cm -3 pc (the DM predicted by the Taylor & 
Cordes 1993 electron density/distance model is ~ 70 cm -3 pc). 
Each of the 204 resulting time series was then searched for peri- 
odic signals over a range of duty cycles with an FFT-based code 
(described in detail by Lorimer et al. 2000) identifying signif- 


icant features in the fundamental amplitude spectrum as well 
as in spectra with 2, 4, 8, and 16 harmonics folded in. The 
1375 MHz data were heavily corrupted by RFI — in particular 
the 182nd harmonic of a radar signal repeating every — 12 s co- 
incided with the expected pulsar rotation frequency. Therefore 
we eventually reanalyzed these data while blanking the radar 
frequency and its first 150 harmonics in the amplitude spectra. 

The 820 MHz data were clean in the spectral range of interest 
and a signal was unambiguously detected at the expected pe- 
riod and clearly dispersed, peaking at DM « 143 cm -3 pc with 
a signal-to-noise ratio S/N = 10.9, as shown in Figure 1. After 
excision of the radar-related RFI, the 1375 MHz data yielded 
a signal at the same period within the uncertainty peaking at 
DM & 140cm" 3 pc with S/N = 8.8. 



Fig. 1 . — Search signal-to-noise ratio as a function of trial dispersion mea- 
sure (DM) for the candidate with period P % 65.69 ms in the 820 MHz data set. 
Each curve corresponds to the addition of the number of harmonics indicated. 

We then measured pulse times-of-arriva! (TOAs) referenced 
to the observatory atomic time standard from both data sets: 6 
good TOAs were obtained in each of the 800 and 1400 MHz 
bands with individual uncertainty % 0.5 ms. A further 17 TOAs 
of similar quality were obtained on 8 additional days of obser- 
vations, with a total data span of 35 d. We used the TOAs and 
the TEMPO timing software 11 to fit for F, F, and DM, obtaining 
the values listed in Table 2. The table also lists the celestial co- 
ordinates of PSR J0205+6449 obtained by Slane et al. (2002), 
which were held fixed in our timing fit. Our values of F and 
F determined from a phase-connected timing solution spanning 
35 d are much more precise than, but agree with, the equiva- 
lent values obtained by Murray et al. (2002) based on only 2 
periods measured 1 .2 and 4.4 yr before: the barycentric period 
F listed in Table 2 differs by 15 ± 18 ns from that extrapolated 
to the epoch of our observations using the F and P determined 
by Murray et al., while our F is consistent with theirs at the 
2cr level. Such a close match of spin parameters over an inter- 
val greater than 4yr suggests that no significant period glitch 
occurred during this time span. 

Both search data sets, de-dispersed at DM =141 cm -3 pc and 
folded at the best search periods, are shown in Figure 2. The 
average pulse profile at both frequencies consists of a single 
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narrow pulse of width ^ 0.05P (see Table 2); such a profile is 
expected based on Figure 1, where the search S/N is a maxi- 
mum for 8 harmonics summed. There is no evidence for the 
existence of an interpulse at either frequency with amplitude 
> 15% that of the sharp pulse, in contrast with the X-ray pro- 
files that show an interpulse with amplitude ~ 30% that of the 
main pulse (Murray et al. 2002). 



Fig. 2. — Pulse profiles of PSR J0205+6449 displayed as a function 
of time ( bottom ) and summed (top). Left : Data at a center frequency of 
v = 1375 MHz presented with 100 phase bins. The pulse full-width at half 
maximum (FWHM) is 2.3 ± 0.3 ms. Right: Data at u - 820 MHz displayed 
with 60 bins. The pulse FWHM = 3.8 ± 0.4ms. All data were de-dispersed 
at DM =141 cm -3 pc and one phase bin in each panel matches the dispersion 
smearing resolution of the respective profile. Two full periods are shown at 
each frequency for clarity, The increasing/decreasing slope of each reduced 
X 2 curve traces pulsation significance as measured against a constant flux (or 
non-pulsed) model. Decreases in reduced x 2 are most likely caused by inter- 
stellar scintillation or RFI. 


Note that the S/N in the 820 MHz data appears to be smaller 
near the beginning of the integration and again towards the end 
(Fig. 2), possibly as a consequence of scintillation. Assuming 
that the search S/Ns in the 2 long integrations reflect the respec- 
tive long-term time-averaged flux densities in both frequency 
bands, we have integrated the area under the pulse profiles and 
measured their baseline rms to determine the respective flux 
densities. We converted to a Jy scale using the observing pa- 
rameters listed in Table 1 together with the radiometer equation 
(Dicke 1946). In this manner we estimate that S1400 ~ 45 /i Jy 
and Sgoo ~ 130//Jy, with ~ 25% accuracy. The resulting spec- 
tral index in the 800-1400MHz range is a ~ -2.1 ± 0.6, and 
the radio luminosity is L] 40 o ~ 0.5 (d/3.2 kpc) 2 mJy kpc 2 . These 
values will be much improved by the averaging of several cali- 
brated observations. 

The 947 yr-old Crab pulsar was detected initially through its 
individual “giant” pulses rather than through its time-averaged 
pulsed emission (Staelin & Reifenstein 1968). As the emission 
of giant pulses may be related to the magnetic field strength 
at the pulsar light cylinder, Z?/ r (e.g., Cognard et al. 1996), and 
PSR J0205+6449 has the 1 0th largest value of B !r among 
known pulsars, we searched for aperiodic, dispersed pulses 
from PSR J0205+6449. No excess of single pulses was found 


at the DM of the pulsar. Inferring a quantitative limit is com- 
plicated by significant levels of RFI, and this search will be 
reported in detail elsewhere. 


4. DISCUSSION 

We have discovered radio pulsations from the direction of 
PSR 10205+6449 in SNR 3C58 in an unbiased search of 2 
independent data sets acquired at center frequencies of 820 
and 1375 MHz. The pulsations are clearly dispersed and have 
a precisely determined period matching that observed from 
PSR J0205+6449 at X-ray wavelengths (Murray etal. 2002). 
We have therefore discovered radio pulsations that are unam- 
biguously from PSR J0205+6449 which supplants the Crab pul- 
sar as the youngest neutron star with known radio emission. 

The DM = 141 cm -3 pc of PSR J0205+6449, located at 
(l,b) = (130.°72,3°08), is approximately twice that predicted 
by the Taylor & Cordes (1993) model for the usually ac- 
cepted distance to 3C58 (3.2 kpc, determined from H I absorp- 
tion measurements 12 ). This implies an average free electron 
density to 3C58 of n e « 0.044 cm" 3 . The neutral column den- 
sity to the SNR is 3 x 10 21 cm -2 (Helfand, Becker, & White 
1995), implying an ionized fraction of roughly 15%, which is 
somewhat higher than might be expected towards the Galactic 
plane. However, the SNR lies on the eastern edge of the large 
W3/W4/W5 complex of Hu regions (e.g., Reynolds, Sterling, 
& Haffner 2001), and in fact may be on the edge of a cavity 
being ionized by the star-forming regions (R. J. Reynolds and 
G. Madsen, private communication). It is not known if the pro- 
genitor of SN 1181 belonged to, or was a runaway star from, 
this region of star formation. In any case, the relative locations 
of 3C58 and W3/W4/W5 may account for the apparent excess 
of ionized material towards the pulsar. 

We note that the reported detection of PSR J0205+6449 
with DM = 24 cm -3 pc at a frequency of 1 10 MHz (Malofeev, 
Malov, & Glushak 2001) is incompatible with the results we 
present: the dispersion smearing timescale for a signal with 
DM= 141 cm -3 pc within one of the 20 kHz-wide spectral chan- 
nels used by Malofeev et al. (2001) at 1 10 MHz is 17 ms. De- 
dispersing and summing their 32 individual such channels with 
DM = 24cm -3 pc would result in an overall smearing of the 
pulse far in excess of P = 65 ms, rendering pulsations unob- 
servable from PSR J0205+6449. 

The main X-ray pulse from PSR J0205+6449 has es- 
sentially the same sharp structure and width (^ 0.04P; 
Murray et al. 2002) as the radio pulse (Fig. 2 and Table 2). This 
morphology suggests that the X-ray emission is likely magne- 
tospheric in nature, but absolute alignment of radio and X-ray 
profiles will provide further clues for understanding the emis- 
sion mechanism(s). Considering the obvious existence of mag- 
netospheric emission, we might expect this pulsar to be a 7-ray 
source. The flux sensitivity of the EGRET instrument to emis- 
sion from PSR J0205+6449, using a photon index T = 2 typi- 
cal of 7-ray-emitting pulsars, is F 1 ~ 5 x 10 32 erg s -1 kpc -2 (> 
100 MeV; Hartman et al. 1999), corresponding to a luminosity 
limit Ly E F^ 2 < 5 x 10 33 ergs -1 . With spin-down luminos- 
ity E - 2.7 x 10 37 ergs -1 , its efficiency at converting rotational 
power into beamed high-energy 7-rays is ?/ = Ly /E < 2 x 10~ 4 
This compares to 77 ~ 10 -4 for the Crab pulsar, with the largest 
value of E known in the Galaxy, and larger efficiencies for 


12 Depending on data set and Galactic rotation curve used, the inferred distance ranges from 2.6 kpc to “just beyond” 3.2 kpc (see Roberts et al. 1993, and references 
therein). We adopt d = 3.2 kpc in this Letter. 
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(generally older) pulsars with smaller values of E such as Vela 
(tj ~ 6 x 10" 4 ) and Geminga (77 ~ 0.02). If PSR J0205+6449 
follows the approximate rj a relation that seems to ap- 
ply to known 7-ray-emitting pulsars (Thompson et al. 1999), 
its flux may be just below the detection threshold of EGRET 
and should be detected by next-generation observatories such 
as GLAST, 

The discovery of extremely weak radio pulsations from 
PSR J0205+6449 highlights the vexing issue of the true lu- 
minosity distribution of (particularly young) pulsars. The past 
year has witnessed the discovery of 3 very young (r r < lOkyr) 
radio-emitting pulsars with 0.5 < L1400 £ 3mJy kpcC 1-2 or- 
ders of magnitude less luminous than any young pulsars previ- 
ously known (Halpem et al. 2001; Camilo et al. 2002, this Let- 
ter). It now seems abundantly clear that we have yet to probe 
the depths of the luminosity function of young pulsars. The 
state-of-the-art for finding young radio pulsars (~ 10 hr integra- 
tions at telescopes with gain > 1 K Jy _1 and bandwidths 100- 
300 MHz at v ~ 1400 MHz) yields limits usually no better than 
L]4 oq ~ 1 mJy kpc 2 . Nevertheless, many worthy targets remain 
to be searched at this level. Improving the existing limits helps 
to constrain better the “beaming fraction” for young pulsars 
(e.g.. Frail & Moffett 1993, based in part on the non-detection 
of the pulsar in 3C58 with a limit of 0.15 mJy at 1400MHz, 
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estimated a beaming fraction of 0.6 which needs to be revised 
in light of recent discoveries), and is important for building an 
accurate census of Galactic neutron stars. In this regard, the 
detection of PSR J0205+6449, an 820 yr-old pulsar, at a fre- 
quency of 820 MHz — only the 2nd pulsar ever discovered in 
the 800 MHz band (D’Amico et al. 1988) — suggests that the 
availability of the magnificent GBT provides a new region of 
phase space with significant, and as yet untapped, discovery 
potential. 
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Table 1 

Summary of Observations and Analyses at GBT 
Observing frequency (MHz) 1375 820 


Telescope gain (K Jy” 1 ) 2.0 2.0 

System temperature on cold sky (K) ... 20 25 

Galactic synchrotron temperature (K) . 2 9 

SNR 3C58 temperature (K) 64 67 

Spectral channel width (MHz) 1.4 0.5 

Sample time (^.s) 50 72 

Integration time (hr) 8.75 8.25 

Effective time resolution of search (ms) 0.9 0.864 

Effective integration time of search (hr) 8.4 8.1 


Table 2 

Parameters of PSR J0205+6449 


Parameter Value 


R.A. (J2000) a . 02 h 05 m 37?92 

Decl. (J2000) a . +64°49'42"8 

Period, P (ms) 65.68638162(2) 

Period derivative, P 1.9393(4) x 10 -13 

Epoch (MJD [TDB]) 52345.0 

Dispersion measure, DM (cm -3 pc) . . . 140.7(3) 

Pulse FWHM at 1 375 MHz (ms) 2.3 ± 0.3 

Pulse FWHM at 820 MHz (ms) 3.8 ± 0.4 

Flux density at 1375 MHz, Suoo (a*J y ) ~ 45 

Flux density at 820 MHz, Sgoo (juJy) . . ~ 130 

Distance of SNR 3C58, d (kpc) ^ 3.2 

Derived parameters: 

Radio luminosity, L^oo (mJy kpc 2 ) . . ~ 0.5 

Spectral index, a ~ -2.1 ±0.6 

Mean free electron density, n e (cm -3 ) & 0.044 


Note. — Numbers in parentheses represent 1 a uncertainties 
in the least-significant digits of fitted parameters. 

“Position known with ^ 0"5 accuracy from Chandra data 
(Slane et al. 2002). 





